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I. Introduction 

If glassy materials are quenched from above their glass transition temperature deep into 
the glassy state, they usually do not reach equilibrium on experimentally accessible time 
scales. The dynamics in this out-of-equilibrium situation, the so-called aging phenomena, 
have been investigated for a long time, for reviews see e.g. refs.p]|2|. If the system reaches 
equilibrium, i.e. if the quench is not performed at a too low temperature, the aging behavior 
often can be well described in terms of the Tool-Narayanaswami-Moynihan model [3]. This 
model takes into account the non-linearity and the time-dependence of the relaxation time 
via the introduction of a reduced time and a Active temperature. 

In the recent past, many theoretical investigations of the aging behavior of glassy sys- 
tems have been undertaken and some rigorous results regarding both, the detailed time- 
dependences of two-time correlation functions (CFs) and the violations of the fluctuation- 
dissipation theorem (FDT) have been derived for several models[4j. These violations can in 
favourable cases be used for the definition of a so-called effective temperature Also some 
computer simulations on model supercooled liquids have been performed and an effective 
temperature could be definedjHJ IZj- Apart from theoretical considerations a number of 
experimental investigations of the violations of the FDT have been performed in the recent 
past, in particular on glassforming liquids 8 , polymersjHj, colloidal glassesfTHj as well as 
spin glasses^!]. In all these systems, strong violations of the FDT have been observed and 
in some cases an effective temperature as a function of the time elapsed after the initial 
quench into the glassy phase has been determined. For glass forming liquids, the effective 
temperature relaxes to the bath temperature for not too deep quenches [HJ Ej- 
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Without referring to an effective temperature explicitly, it has been found experimen- 
tally that the aging behavior of volume and enthalpy in a glassy polymer may not be the 
same[I2]. A detailed study of the aging properties of polystyrene has been performed by 
Thurau and Ediger[T3]. who clearly showed that probe rotation and probe translation show 
different behaviors in this example. It was found that after a temperature jump into the 
glassy phase, the relaxation towards equilibrium of probe rotation takes longer than the 
corresponding re-equilibration of probe translation, if the final temperature of the quench 
is low enough. In particular, the temperature dependence of the re-equilibration times for 
rotations was found to be stronger than that of translational motions. The results could 
consistently be interpreted in terms of spatially heterogeneous dynamics, which is now 
well established to be at the origin of many of the peculiar features of the relaxation in 
glass forming liquids, for reviews see refs. [TH ITS] . Thurau and Ediger did not attempt to 
quantitatively discuss the relation between the re-equilibration times and the time scale of 
the primary relaxation at the final temperature of the quench. Such an analysis has been 
performed in detail for a number of glass forming liquids by dielectric spectroscopy JH] and 
it was shown unequivocally that the aging behavior in these examples is determined by the 
structural a-relaxation. 

In the present paper, I will discuss the aging properties of glass forming liquids within 
a simple free-energy model that has been used previously to describe the heterogeneous 
relaxation in equilibrium |17l H%1 119j. The idea of this model is simple. It is assumed that 
the primary relaxation in a supercooled liquid is thermally activated and proceeds via 
transitions among a large number of metastable states, characterized by their free-energy. 
This idea is in accord with various proposed scenarios of glassy relaxation, such as the one 
of a random first order transition [20] . In order to describe the rotational or translational 
motion of tagged particles, it is assumed that each transition is accomponied by a particle 
rearrangement, giving rise to changes in both, the orientation and the position of the par- 
ticle considered. These assumptions along with a mean jump distance and a mean rotation 
angle are enough to compute relevant CFs, such as the intermediate scattering function 
or rotational CFs. Due to the different ways various CFs average over the distribution of 
free-energies a number of hitherto unresolved puzzling features of the a-relaxation could be 
explained in a simple wav|18j. This includes the breakdown of the Stokes-Einstein relation 
and the differences in the stretching observed in different experiments probing the reori- 
entational or translational dynamics. Here, I will concentrate on a qualitative description 
of aging experiments and show that the results from model calculations capture the most 
prominent features observed experimentally. 

The paper is organized as follows. In the next section, the free-energy model is briefly 
reviewed and it is described how to calculate CFs in an aging system. In Section III 
the aging behavior of the intermediate scattering function and of the rotational CFs is 
discussed. For comparison, in Section IV the aging-behavior of rotational CFs is discussed 
in the framework of models for environmental fluctuations |21j or exchange models [22 . 
Finally, the conclusions are presented in Section V. 
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II. Aging behavior in the free-energy landscape model 



The aging behavior of supercooled liquids depends on the final temperature in a tempera- 
ture jump experiment. This is because if the quench is performed to a temperature not too 
far below the glass transition temperature, the system can reach equilibrium on an acces- 
sible time scale since the glass transition temperature is merely a convenient way of classi- 
fying the time scale of the primary relaxation |23j. For deep quenches the re-equilibration 
usually cannot be monitored experimentally because the relaxation time is extremely long. 
This is one reason why sometimes models exhibiting a true transition into a glassy phase 
are considered when discussing the aging behavior. As will become clear in the present 
Section, in the free-energy landscape model the system always reaches equilibrium for long 
times and the re-equilibration time scale is determined by the a-relaxation time at the final 
temperature of the quench. 

In the following, I will utilize the free-energy landscape model for the caculation of 
translational and the rotational two-time CFsjTH| after a temperature jump. In case of the 
reorient at ional motion, the CF is obtained by correlating the orient at ion- dependent inter- 
actions relevant in the experiment considered at two times. The transformation properties 
of these interactions determine the so-called rotational CF gL(t,t w ), 

g L (t,t w ) = (P L (cos(0(*))P L (cos(0(O)> (1) 

Here, Pl{x) denotes the Legendre polynomial of rank L. For instance, in dielectric spec- 
troscopy g\(t, t w ) and in NMR or depolarized light scattering g 2 (t, t w ) are measured. Here, 
t w denotes the waiting time, i.e. the time that has elapsed after the initial preparation 
of the system before the beginning of the experiment. In all later calculations, the initial 
preparation corresponds to a quench from a high temperature to the working temperature. 

If one is interested in translational motions of tagged particles, one naturally considers 
the incoherent intermediate scattering functionjM: 

S q (t,t w ) = (e^We-M*.)) (2) 

where r(i) denotes the position of the particle and q is the wave- vector. Additionally, I 
have already assumed that the system is isotropic, q= |q|. 

In the free-energy model, a free-energy e is associated with each of the exponentially 
large number of basins or metastable states of the system jH The a-relaxation is 

modeled using a master equation j2H] for the conditional probability to find the system in 
a metastable state characterized by e at time t, given that one had eo at to: 

G(e, t\e , t ) = - J dtK(e\e)G(e, t\e , t Q ) + J deK(e\e')G(e, t\e , t Q ) (3) 

Here, K(e'\e) denotes the rate for a e — > e'-transition, which is assumed to be time- 
independent for simplicity. Consequently, G(e, t|e , to) depends only on the difference of the 
times involved, G(e, t\e , t ) = G(e, t — t 1 e ) . In ref.JTB], we considered different choices for 
the transition rates ft(e'|e) with only quantitative differences in the results. Therefore, in 
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the present calculations, I will solely consider the following simple choice for the transition 
rates, denoted as globally connected model in ref. 18 . It is assumed that every state can 
be reached by a transition and that the destination state is chosen at random according 
to a prescribed density of states, denoted by 77(e). The transition out of a given state is 
assumed to be thermally activated with a rate given by n(e) = K 00 e^ e ) where a common 
activation energy has been absorbed in the prefactor and (3 = (fe^T) -1 . This means, all 
calculations will be restricted to the simple case of transition rates 

<e'\e) = /wrtOe* (4) 

Furthermore, the density of states will be chosen to be Gaussian with zero mean and 
variance a throughout, 77(e) = i- e~ e ^ 2a2 \ The same transition rates have also been 
used by Dyre in his 'energy master equation' |27| . With this choice, the system reaches 
equilibrium for long times and the equilibrium probabilities are given by: 

p eq {e) = lim G(e,t\e ,t ) = -^^ e - {e - e ~ )/(2,j2) with e = -(3a 2 (5) 

The distribution of equilibrium probabilities thus has a temperature-independent width 
and a mean that scales like the inverse temperature. In the inherent structure approach to 
the classification of the potential energy landscapes of supercooled liquids f.25] it is usually 
found that the inherent structure energies are distributed according to a Gaussian [28 . 
However, the transitions among different inherent structures are more complicated than 
assumed in eq.flU), see e.g. ref.[29j. 

The main idea of the present model consists in the coupling of the reorient at ional and 
translational motion of a tagged particle to the a-relaxation as follows. The transitions 
among the metastable states are the only possible way for particle rearrangements (apart 
from vibrational motion within a state). Therefore, it is assumed that the particle con- 
sidered changes its position as well as its orientation by a certain amount whenever a 
transition takes place. In order to model the dynamics of either translations or rotations, 
the composite Markov process {e(t),a(t)} with a(t) denoting the orientation Q(t) or the 
position r(i) is consideredjd [THj. To this end, one has to deal with the master equa- 
tion for the combined conditional probabilities G(e, a, t\eo, ao, to) and the transition rates 
W(e',a'\e,a). For the latter, the simplest possible choice will be considered which means 
that only reorientations with a mean angle 5Q and only translations with a mean jump 
distance 5R are allowed. 

The solution of the corresponding master equation is outlined in Appendix A for conve- 
nience of the reader, for more details see ref.^Sj- Here, it suffices to note that it is assumed 
that initially there are no correlations between the populations p(e, 0) and p(a, 0) and the 
latter are given by the equilibrium populations p eq (a), cf. eq. fTOjl . The two-time CFs 
defined in eqns.((T} and (j2j) are given by: 

g^t, t w ) = JdeJ de'G L (e, t\e')p(e', t w ) 

S Q (t,t w )= JdeJde'G Q (e,t\e')p(e',t w ) (6) 
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Here, the Green's functions G^(e, t\e') and Gq(€, t\e') are solutions of eq. (jA.5|) and Q = 
(q-SR). For long waiting times, it is evident from eq.(0) that 

p(e,t w ) ->p eq (e) for t w -> oo (7) 

because of p(e, t w ) = Jde'G(e, t w \e')p(e', 0) and Jde'p(e',0) = 1. The CFs then reduce to 
those in equilibrium, 

gl\r) = JdeJde'G L {e,r\e')f\e') and S%(r) = J dej de'G Q (e,r\e')p e \e') (8) 

the former of which has been discussed in detail in ref.|18j. 

In order to describe an aging experiment, one has to fix the initial populations. Ex- 
perimentally, one usually performs a temperature jump from a 'high' temperature above 
the glass transition temperature T g to below T g . In the present paper, I will consider only 
quenches from /3 = (T = oo) to the working temperature, which means: 

p(e, 0) = V {e) = -^e~^ (9) 



Due to the fact that equilibrium is reached in the long run, cf. eq.(J7J), all aging effects 
are of a transient nature. Note that (3 = in the present context only means that the 
transition rates in eq.(jU are no longer thermally activated, i.e. one has K^ =0 ^(e'\e) = 
K ooV( e ') independent of the initial state. 

A first estimate of the effects that are to be expected can be obtained from the temporal 
evolution of the populations, p(e,t w ). In Fig. la, I have plotted p(e,t w ) as a function of e 
for two final temperatures, T = 0.2a (upper panel) and T = OAa (lower panel), for various 
waiting times t w . For t w — 0, one has a Gaussian centered at e = 0, cf. eq.© and in the 
long-time limit, the center has moved to e = —f3a 2 . It is, however, obvious from the figure, 
that with increasing t w the distribution first narrows and then becomes broader again. 
Furthermore, the form of the distribution is not exactly Gaussian for intermediate times. 
A similar behavior of the distribution of inherent energies during aging has been found in 
a recent simulation of a binary Lennard- Jones system 30J. In a next step, I determined 
the various moments of the distributions, the mean value (e(t w )), the width a(t w ) and the 
skewness ji(t w ). The results are shown in Fig. lb. The change in sign of 71 (t w ) reflects the 
fact that the distributions are skewed to the high-energy side for small t w and to the low- 
energy side for long t w . The deviations from a Gaussian are seen to be more pronounced for 
the lower final temperature. The fact that the distribution of p(e, t w ) becomes narrower for 
intermediate t w already gives rise to the expectation that the degree of non-exponentiality 
of the two-time CFs given in eq.(JH} will change with the waiting time. This means that 
one expects the stretching parameter [3k of Kohlrausch fits (oc exp [-{t/T K Y K }) first to 
increase and then to decrease again. 

In the following, the evolution of the intermediate scattering function and the rotational 
CFs after a quench from high temperature to several low temperatures will be discussed. 
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III. Aging monitored by two-time correlation functions 



The equilibrium properties of the rotational CFs g^{r) f° r different values of L have 
been discussed in detail in ref. [IS] along with the apparent translational enhancement. 
Therefore, I will start with a brief discussion of the Q-dependence of the intermediate 
scattering function. In Fig. 2a, Sq{t) is plotted versus t/tq for various values of Q and 
T = 0.3a. Here, tq is the relaxation time obtained from a fit to a stretched exponential 
function exp (—(t/tq)^ q ). It is evident that only for small Q one is in the hydrodynamic 
limit and the degree of non-exponentiality is a strong function of Q. In the free-energy 
landscape model, the reason for this behavior is explained by the fact that for small Q 
many jumps have to be performed in order for the CF to decay. Therefore, Sq(t) averages 
over the density of states quite effectively and the decay is very slow and exponential. For 
large values of Q, any single jump decorrelates the intermediate scattering function and 
therefore there is no averaging. Consequently, the decay is governed by a broad distribution 
of effective relaxation rates. The situation thus is similar to the one for the rotational CFs 
for different values of L, c.f. the discussion on this point in ref. ^H]- For rotational CFs, 
however, the mean rotation angle enters separately, whereas due to Q — q-5R the jump 
distance 5R only multiplies the scattering vector. 

In Fig. 2b the relaxation times tq and the stretching parameter (3q from fits of Sq(t) 
to a stretched exponential function, exp (— (t/tq)^ q ), are shown. For values of Q smaller 
than roughly 10~ 2 , the decay is exponential and one can define a diffusion coefficient 
D t = (5R) 2 (k}/6, where (k) is the average relaxation rate, (k) = JdeK(e)p eq (e). For larger 
Q, @q decreases and the relaxation time deviates from the Q _2 -scaling. For Q > 10, neither 
tq nor (3q change any further, because here one is in the limit of the 'jump correlation 
function' [T7j . 

In Fig. 2c, I have plotted the fitted relaxation times scaled to those at T = 0.9a versus 
inverse temperature. The reason for the different behavior again lies in the different av- 
eraging over the transition rates for different Q. The maximum discrepancy between the 
relaxation times, however, does not exceed a factor of roughly two (three, if the average 
times are considered). The values of Q were chosen, because if 5R is chosen between twenty 
and fifty percent of the van der Waals radius of a tetracene molecule, they are comparable 
to the scattering vector used in the mentioned experiment by Thurau and EdigerfTB*!. 

Next, the CFs Sg(t w + T,t w ) and + t, £„,) are calculated according to eq.©, where 
now t = t w + t is used as the time- variable. Results for different waiting times t w are shown 
in Fig. 3a for T = 0.3a. In case of the rotational CF, I used L — 1, corresponding to 
dielectric relaxation. The results for other values of L are quite similar, apart from minor 
quantitative differences that are relevant at equilibrium, cf. ref.^B]- ^ n a ^ calculations, 
a mean rotation angle of = 10° is used. This is on the order of magnitude as found 
experimentally for supercooled liquids |3Tj. For long t w , both CFs are independent of t w 
and coincide with those in equilibrium. It is evident, that the relaxation time as a function 
of t w changes much more for the rotational CF than for the intermediate scattering function. 
Also this finding is explained by the different inherent averaging over the density of states. 
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As already pointed out above, for small Q, many e — > e'-transitions have to take place in 
order for Sq to decay. The same transitions are responsible for the approach of equilibrium. 
Therefore, the t w dependence of tq is less pronounced than the corresponding one for n, 
the relaxation time decribing g\. This is because in the latter case fewer e — ► e'-transitions 
are required for g% to decay and the system exhibits less re-equilibration during the decay. 

In Fig. 3b, the results of Kohlrausch fits exp(— (t/r z (t M ,)) /3z ^)), where z denotes either 
L— 1 or Q, to the CFs as shown in Fig.3a are collected. Two features are evident immedi- 
ately. The values of T\(t w ) change by about 1.5 decades from t w = to t w = oo, wheras this 
change is much smaller for Tq(t w ) (0.7 decades for Q = 0.1 and 0.07 decades fo (5 = 0.01). 
The stretching parameters show an increase for small t w and then decrease as equilibrium 
is reached. This 'hump' for intermediate t w has its origin in the narrowing of the distri- 
bution of the populations p(e,t w ), cf. Fig.l. Thurau and Ediger[T3| found that r TOt (t w ) 
changes by a larger amount as a function of t w than T tTanB (t w ) does. In their interpretation 
of the experimental results, the value of the stretching parameter for the rotational CF 
was assumed to also show a slight increase followed by a decrease as a function of t w . Both 
effects are much more pronounced in the present model calculations, which might have its 
origin partly in the high initial temperature used here. It should be mentioned, that the 
ratio T\{t w = oo)/ri(t w = 0) strongly depends on the final temperature of the quench, for 
instance this ratio is roughly 15 and 6 for T = 0.32<r and T = 0.35<r, respectively. As al- 
ready mentioned, it will not be attempted to provide a quantitative description of existing 
experimental data. The overall qualitative features are, however, very similar to what is 
observed experimentally. 

Another important feature of the experiment, which supports the view of a spatially 
heterogeneous aging scenario is the fact that probe rotation needs longer to reach equilib- 
rium than probe translation, at least for low final temperatures|13j. In Fig. 3c, T\(t w ) and 
T Q{tw) are plotted in a scaled way. It is evident that the intermediate scattering function 
reaches equilibrium much faster than the rotational CF. In order to have a measure of 
the re-equilibration time, to be denoted as r eq , I have determined the (1 — l/e)-points of 
the curves, cf. Fig. 3c. These are plotted versus inverse temperature in Fig. 3d, where I 
scaled them to their value at T = 0.5a. It is seen that the temperature-dependence of the 
re-equilibration times is somewhat different for translation and rotation. For lower temper- 
atures, T eiJ (rot.) is longer than r eq (trans.) in qualitative agreement with the experimental 
results. The explanation of this finding again lies in the fact that the two CFs average 
in different ways over the density of states. Therefore, Sq{t w + r, t w ) reaches equilibrium 
faster than gi(t w + r,t w ). It should be pointed out that in the framework of the present 
model it is the same difference in the averaging over the density of states that gives rise to 
the apparent translational enhancement |18j . Furthermore, the re-equilibration of course is 
determined by the a-relaxation in the free-energy model as there is no other time scale. 
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IV. Aging in environmental fluctuation models 



In the free-energy model the life-time of the dynamic heterogeneities is on the same time 
scale as the a-relaxation time. This is because there is no other time scale in the model. In 
terms of Heuer's more general concept of a rate memory parameter Q (not to be confused 
with a scattering vector) the free-energy model intrinsically corresponds to Q = 1 [32] • It has 
to be noted, however, that the free-energy model in the simple form presented here, cannot 
describe the finding of an exchange time much longer than the a-relaxation time observed in 
a photobleaching experiment near T g in ortho-terphenvl [3lH I15j. In this section, therefore, 
I will consider so-called environmental fluctuation models [2TJ [22] that allow for an extra 
time scale for exchange or arbitrary rate memory. In such models it is assumed that the 
molecular dynamics are faster in some part of the sample and slower in another. When 
applied to supercooled liquids the dynamic heterogeneities are related to a distribution 
of the correlation times and their life-time is limited by exchange among the regions of 
different mobility. Therefore, one would interpret the dynamic exchange as the structural 
relaxation of the system. 

Such models can also be used for the calculation of the aging properties of the relevant 
CFs. Here, I will concentrate on the rotational CF gi(t,t w ) and consider the following 
variant of an exchange model. It is assumed, that a reorientation rate T(e) corresponds to 
each value of the variable e. Therefore, the variable e in the present context characterizes a 
given environment and the value of the corresponding reorientation rate T(e). These rates 
for simplicity are chosen according to 



The transition rates n(e\e') for e' — > e-transitions are again chosen according to eq.(jU). As 
mentioned above, it are these transition rates that are responsible for the exchange among 
the various reorientation rates in the sense that with every e' — ■> e-transition a change from 
T(e') to T(e) is accomponied. According to eq.©, one has a Gaussian distribution of 
e- values and therefore a broad distribution of reorientation rates T(e). This 'bare' distribu- 
tion is narrowed somewhat due to exchange, if the latter takes place on a similar time scale 
as the reorient at ional dynamics. Note that the choice of T(e) and n{e\e') are of a purely 
phenomenological nature. However, the details of the form of the distributions of reori- 
entation rates and of exchange rates are of minor importance for the following qualitative 
discussion. 

The solution of the master equation for such exchange models is outlined in Appendix 
B where again it is assumed that reorientations proceed via rotational jumps with a mean 
jump angle = 10°. In the following, I will restrict myself to the assumption that each 
'exchange-transition' is accomponied by a randomization of the molecular orientation. In 
this case, the solution of eq. flB.2j) is trivial (because of cl = 6l,o) and the corresponding 
rotational CF is given by gL(t,t w ) = Jdee~ ( - FL ^ +K ^^ t ~ tw ' > p(e,t w ) with the reorientation 
rates Ti(e) = [1 — Pl(cosO)] T(e), cf. eq. flB.3j) . The aging behavior of gi(t,t w ) is reflected 
solely in the populations p(e,t w ) which evolve in time according to the master equation 




(10) 
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(J3J). These populations decay with effective reorientation rates (T L (e) + n(e)) , reflecting the 
fact that both, a reorientation as well as an 'exchange-transition' depletes the population 
p(e,t w ). In order to allow for different time scales of reorientations and exchange, I will 
use the following simple relation between the attempt frequencies: 

^ = x with r'^ = [1-cosGjr^ (11) 

Here, x is a free parameter and the scaling to (1 — cosG) assures that the reorientational 
geometry has no further impact on the relative time scales. For x — 1, the time scale of 
molecular reorientation coincides with the time scale of dynamic exchange, whereas for 
larger x, the average exchange rates are smaller than the reorientation rates. Note that 
the assumption of a randomization of the molecular orientation means that also exchange 
gives rise to a decay of the rotational CF, which has some influence on its form for x ~ 1 
but not for large x (slow exchange). 

The rotational CFs gi(t,t w ) for x—1 behave very similar to those obtained in the free- 
energy landscape model. For large x, however, the re-equilibration takes very long. This is a 
consequence of the fact that the aging behavior is determined by the time scale of exchange 
in this model. The e — > e'-transitions are responsible for both, the re-equilibration of the 
system and the exchange among the reorientation rates. For large x, the re-equilibration is 
still determined by the exchange time but the time scale for reorientations is faster. This 
behavior is summarized in Fig.4, where I have plotted the rotational correlation times T\ 
as a function of the waiting time for T = 0.4a and = 10°. The re-equilibration for x = 10 3 
is much longer than for x—1. This behavior scales exactly with x and it shows that in 
this model the re-equilibration is determined by the dynamic exchange and not by the 
reorientations. 

These model calculations show that within such an approach the assumption of an 
exchange time much longer than the a-relaxation, as obtained in the mentioned optical 
experiments (221 near T g , cannot be reproduced by assuming that dynamic exchange 
is governed by a structural relaxation which also determines the aging behavior of the 
system. This is because the re-equilibration time scale has been found to be determined 
by the a-relaxation^E]- Therefore, one would have to include still another time scale, an 
'aging-time', in such exchange models. This 'aging-time', however, would have to coincide 
with the a-relaxation time. 

As already noted above, also the free-energy model in the simple form presented here 
cannot describe both features. It is not clear at present, whether the reason for this lies in 
the oversimplified assumptions regarding the geometry of molecular reorientations and if 
a more realistic description using a distribution of jump- angles jM] would allow to obtain 
longer exchange times. 
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V. Conclusions 



In the present paper I have examined the out-of equilibrium behavior of the two-time 
CFs for translational and rotational degrees of freedom within the framework of a free- 
energy model for relaxation in glassforming liquids. The system has been prepared in 
an off-equilibrium situation by performing a quench from a high temperature to a low 
'working temperature'. Instead of trying to obtain quantitative agreement with existing 
experimental data the focus has been on a qualitative description of the physics inherent 
in the model. Therefore, all model calculations were performed using the simplest choices 
possible for the various parameters involved. For instance, the transition rates were chosen 
according to a simple-minded rule, cf. eq.(JIJ), and the density of states was chosen to be 
Gaussian with a temperature-independent width. Both choices would have to be improved 
in order to obtain quantitative agreement with experimental results. In addition, I solely 
considered the generic situation of a quench from an infinitely high temperature. 

When considering the behavior of the intermediate scattering function and of the ro- 
tational CFs after a quench from a high temperature to a rather low temperature it is 
found that the correlation times of both functions increase continuously as a function of 
the waiting time that has elapsed between the quench and the measurement. For long 
waiting times, both correlation times reach their equilibrium values at the final tempera- 
ture. However, the absolute change in correlation time is much larger in case of rotational 
motion than for translations monitored with small scattering vectors. In addition, the 
time needed to re-equilibrate, i.e. to reach equilibrium at the final temperature, also is 
longer for rotations than for translations. When considered as a function of the final 
temperature, the re-equilibration times r eq behave differently For rotations, r eq shows a 
stronger temperature-dependences than for translational motions. These findings along 
with a waiting-time dependent stretching parameter are in qualitative agreement with 
experimental results by Thurau and Ediger[ 1-T. 

In the present model, the re-equilibration is determined solely by the a-relaxation, as 
there exists no other time scale. This feature is in accord with recent dielectric experiments 
on a variety of supercooled liquids[in]. To allow for two distinct time scales I considered 
a so-called exchange model in addition to the free-energy landscape model. Such a model 
naturally provides two distinct time scales, one for on-site relaxation and one for dynamical 
exchange. It is found in the model calculations that r eq coincides with the equilibrium 
'exchange time' if the exchange is associated with the structural relaxation responsible for 
the aging properties. Therefore, if the exchange time is much longer than the correlation 
time, r eq is much longer, too. This means, in order to reproduce r eq on the order of r a in 
such models, yet another time scale has to be introduced. 

In summary, the free-energy model that has been introduced in order to describe the 
primary relaxation of supercooled liquids in equilibrium [HI HHj can also reproduce the 
experimentally observed features in out-of-equlibrium situations, at least on a qualitative 
level. The comparison with the exchange models shows, that aging experiments might also 
yield additional information about the nature of the dynamic heterogeneities in supercooled 
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liquids. In particular, it appears an open question what the exact relation between the 
life-time of these heterogeneities and the re-equilibration time after a quench is. In the free- 
energy model the differences in the time scales, the correlation times, the life-time of the 
dynamic heterogeneities and r eq can unequivocally be attributed to the different averaging 
over the density of states for different dynamical variables. The underlying intrinsic time 
scale is determined solely by the a-relaxation. 
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Appendix A: Solution of the master equation for the 
free-energy landscape model 



In this Appendix I will give the relevant formulae that are needed for the solution of 
the master equation for the composite Markov process {e(t),a(t)}, where a(t) has to be 
identified with the orientation Q or the position r. Because the solution of the master 
equation proceeds in the same way for both cases, thetreatment can be formulated quite 
generally. The transition rates can be written in the form 

W(e', a'\e, a) = K(e'\e)A(a'\a) (A.l) 

Here, I neglected possible dependencies of A(a'\a) on the initial and final states of the 
transition. In the text, only the special case A(a'\a) = S(a' — (a + 5a)) is considered 
explicitly. To solve the corresponding master equation for this case one defines the matrix 
of the eigenvectors of A(a\a'), U(a,p) corresponding to the eigenvalue A(p), i.e. A(p) = 
Jda Jda'U~ 1 (a' ,p)A(a'\a)U(a,p). Next, the conditional probability is expanded in terms 
of these eigenvectors, 

G(e, a,t\e ,a ) = J dpU(a,p)G p (e,t\e )U' 1 (a ,p) (A.2) 

The Green's functions G p (e,t\eo) are then found from the solution of 

G p (e, t\e ) = -n(e)G p (e, t\e ) + A(p) J dd K{e\e')G p {e' , t\e ) (A.3) 

where the sum rule «(e):= Jde'n(e'\e) was used for the diagonal element. 
The initial populations are chosen as described in the text, i.e.: 

p(e,a,t = 0)=p e »p(e,t = 0) (A.4) 

If one considers isotropic rotational motions with a fixed mean jump angle SQ = Q, one 
has to identify a with Q. In this case one has p eq {a) = U (a, p) = J^^Dj£l(Q), A(p) = 
Pl(cos (G)) and the integration over p is now a sum over L. In case of translational jumps 
onto all postions of a sphere with a radius 5R (jump length), the corresponding substitutions 
in the general formulae are p eq (a) = V~ x with V denoting the volume. Furthermore, one 
has U(a,p) = ^e* qr and A(p) = j (q6R), where jo(x) denotes the Bessel function of zeroth 
order. 

In eq. (|A.3)1 one has to identify p = L for rotations and p = Q with Q = (q- 5R) for 
translations and thus eq. (|A.3J) explicitly reads: 

Rotation : G L (e, t\e ) = -n(e)G L (e, t\e ) + P L (cos (Q)) J de' n(e\e')G L {e' , t\e ) 

Translation : G Q (e,t\e ) = -K(e)G Q (e,t\e ) + j (Q) J de' K(e\e)G Q (e' , t\e ) (A.5) 

Note that due to Pq(x) = jo(0) = 1, there is one eigenvalue, p — 0, A(0) = 1. Therefore, 
for p—0, from eqns. (|A.3[ IA.5j) the original master equation for the transitions in the free- 
energy landscape, eq.©, is obtained, i.e. Go(e, t\eo) = G(e,t\eo). This is in full accord 
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with the idea underlying the model that the variable a is allowed to change only in case 
of a e— >e' transition[18j. Once the G p (e,t\eo) are obtained from the (numerical) solution 
of the equations (|A.5jl . all quantities of interest can be calculated and one finds for the 
normalized correlation function: 

C p (t, t w ) = JdeJ de'G p (e, t - t w \e')p(e', t w ) (A.6) 

where p(e,t w ) = Jde'G(e,t\e')p(e',t = 0), i.e. only p = is relevant here. This is easy 
to understand from the fact that according to the chosen initial conditions, eq. (jA.4|) . the 
variable a was in equilibrium in the beginning, t w = 0. In case of reorient at ional motions 
the normalization eliminates a factor (2L + l) -1 . 

In general, the behavior of the correlation functions is obtained from a numerical so- 
lution of the equations for the Green's functions, eq. ([A.5|) [17l ITS]. For this purpose, the 
transition rates K(e'\e) have to be chosen in a prescribed way, e.g. as done in the text, 
cf. eq.(HJ). It is, however, possible to give some general results valid in specific situations. 
Consider the case of rotational random jumps, for which one has A(p>0) = 0. Similarly, if 
the incoherent scattering function is observed for large wave- vectors, Q^>1 (or q^> (5R) _1 ), 
one has jo(Q) — 0. 

In the limit A(p) — > 0, the general expression for the Green's function, eq. (jA.3)L ap- 
proximately reads as G p (e, t\eo) = —K(e)G p (e,t\eo) and therefore one has G p (e,t\e ) ~ 
5(e — eo)e _K< - e ' )t . Using this expression, one easily finds from eq. (|A.6)) : 

C p (t, t w ) ~ J dee'^^pie, t w ) (A.7) 

which is independent of p and gives the probability that the system has not left the 
metastable state with free-energy e occupied at t = t w in the time interval (t — t w ). 

In the other extreme, namely very small wave-vectors or rotational diffusion, one has 
A(p) = 1 — T] p with t]q = Q 2 /6 and tjl = L(L + l)(6/2) 2 . As has been shown in ref.jlHj, in 
equilibrium in the limit of small wave- vectors the model predicts an exponentially decaying 
intermediate incoherent scattering function, Sq(t) ~ exp (— q 2 Drpr), with an apparent 
diffusion coefficient D?= (5R) 2 (n)/6. In case of the rotational diffusion of molecules one 
accordingly has g e L q (r)~exp (—L(L + l)D R r) with D R = (6/2) 2 (k). 

Appendix B: Environmental fluctuation model 

Such models can be defined by transition rates of the type [22! 

W{e', Q'\e, fi) = W e {Q'\Q)5{e' - e) + K{e'\e)A{Q'\Q) (B.l) 

where now the on-site reorientations are modelled via finite W e (f)|0'). Translational mo- 
tions can be modeled in a similar way. In the following, it is assumed that the reorientations 
proceed via rotational jumps with a mean jump angle as in the free-energy landscape 
model. In addition, it has to be quantified what happens to the molecular orientation 
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in case of e— >e' transition. As in the original model of Beckert and Pfeiferj35], it will be 
assumed that either no change at all, A(f2|fi') = 1, or a random rotation, A(fl|Q') = l/(87r 2 ), 
takes place. Proceeding exactly in the same way as in Appendix A, the Green's functions 
are obtained from: 

G L (e,t\e ) =-\T L {e) + K {e)]G L {e,t\e )+c L J de' K (e\e')G L (e' , t\e ) (B.2) 

where the reorientation rates r^(e) are given by 

r L (e) = [l-P L (cos0)]r(e) (B.3) 

Additionally, I defined c L = 1, if A(Q\£Y) = 1 and c L = 5 Lfi if A(fi|fi') = 1/(8tt 2 ) [221 - As in 
the case of the energy landscape model, the relevant quantities are given by eqns. (|A.6|) if 
again it is assumed that there is no correlation initially. 
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Figure captions 

Fig.l : a: The distribution of populations p(e,t w ) as a function of free-energy e for dif- 
ferent values of the waiting time t w after a quench from T = oo to T = 0.2er (upper 
panel) (t w = 0, 10~ 15 , 10~ 13 , HT 11 , 10~ 9 , 10~ 7 , 10" 5 , 10" 3 , 10" 1 , 10 2 , 10 4 ) and T = 0.4a 
(lower panel) (t w = 0, 10" 5 , 10~ 3 , KT 1 , 10, 10 4 ). 

b: The mean value (upper panel), the width (middle panel) and the skewness (lower 
panel) of the distributions shown in (a). The skewness is defined as ji(t w ) = 
K>3(t w ) / K2(t w ) 3 / 2 where k 2 and n 3 are the time- dependent cumulants. 

Fig. 2 : a: The intermediate scattering function in equilibrium (t w — > oo) for T = 0.3<r 
and different values of the reduced scattering vector Q = q-8R. 

b: Results of fits to a Kohlrausch function (exp (— (£/ t q) /3q )) a ^ T = 0.3a as a func- 
tion of Q. It is evident that Q 2 t q = const, and Pq — 1 holds only for small Q. 
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c: The parameters tq (upper panel) and (3q (lower panel) as a function of tem- 
perature. The values of tq are scaled to their value at T = 0.9a, TQ(T) scaled = 
rg(T)/rg(0.9). 

Fig.3 : a: S Q (t w + T,t w ) and g^ty, + T,t w ) versus r sca i cd = (t/t q ) and r scaled = 
(t/ti), respectively. Both panels show the correlation functions for reduced t w = 
10~ 7 , 10 -1 , 1, 10, 10 6 tq(ti). For the longest waiting time, the system has reached 
equilibrium. 

b: Results of Kohlrausch fits to Sq(t w + r, t w ) and gi(t w + r, t w ) as a function of the 
waiting time. 

c: The time scales TQ(t w ) and Ti(t w ) versus t w are shown in a scaled way such that 
t x (0) =0 and t x (oo) = 1 for x = Q, 1. 

d: The re-equilibration time, i.e. the time needed for the correlation time to reach 
its equilibrium value. The values of r eq are determined as the (1 — l/e)-values of the 
curves shown in part (c). 

Fig. 4 : The rotational correlation time Ti(t w ) versus (t w /ri) for the exchange model 
showing that for x = 10 3 the re-equilibration takes much longer than the rotational 
correlation time. 
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G. Diezemann, Fig. 1a 
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G. Diezemann, Fig. 1b 
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G. Diezemann, Fig. 2a 
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G. Diezemann, Fig. 2b 
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